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I. INTRODUCTION 



It has been widely known that the spectral splitting of the lowest two states of a quantum 
mechanical symmetric double-well potential due to the quantum tunneling is successfully 
calculated by summing up multiple instanton contributions with the dilute gas approxima- 
tion On the other hand, it might not have been duly recognized that a naive application 
of the method to a bit more complicated system in general confronts some new difficulties 
and hardly yields proper results. In this respect, the problem of a symmetric triple- well 
potential described essentially by the following function: 

V{q) = \q\q' - l)^ (1) 

have been recently attracting attention of several research groups. One of the novel features 
of the above potential over the symmetric double-well potential, besides the obvious differ- 
ence in the number of the potential wells, comes from the fact that the harmonic frequency 
of the central well is different from that of the side wells. As a result, there are no degen- 
eracies between the harmonic oscillator spectrum of the central well and the side wells while 
that of the left and right well completely degenerate with each other. Thus, it is difficult 
to expect intuitively how the quantum tunneling effect contributes to the each harmonic 
oscillator spectrum. 

To the best of our knowledge, the multi-instanton calculation technique with the aid of 
the dilute gas approximation was applied to a triple-well potential problem first by Lee et al. 
01 and later independently by Casahorran . Both of their resulting formulas of the lowest 
three eigenvalues are however peculiar and doubtful in the fact that they do not coincide with 
the harmonic oscillator spectrum of the each potential well when the instanton contribution 
due to the quantum tunneling effect is turned off, although both of the authors did not 
discuss the validity of them nor compare the obtained results with the exact eigenvalues. 

A few years ago, we investigated in Ref. a similar problem in a different context, 
namely, dynamical breaking of A/'-fold supersymmetry, by means of the valley method 
which is a generalization of the semi-classical approximation and enables us to calculate 
nonperturbative correction beyond the dilute gas approximation. Our formulas of the spec- 
trum are completely different from the ones previously obtained in Refs. jH |^ and we have 
justified our results by checking consistency with some characteristic features of A/'-fold su- 
persymmetry discussed earlier in Ref. j^. We have found that the contribution from the 
interaction between the instantons plays a crucial role in the calculation in order to yield 
the correct formulas. Hence, we have asserted that the dilute gas approximation would fail 
in the case. 

Recently, however, Alhendi and Lashin reexamined the triple-well potential problem and 
carried out a careful calculation of the multiple instanton contribution with the dilute gas 
approximation 0. Their results are different from ours but look sensible in the sense that 
they reduce to the harmonic oscillator spectrum when the instanton corrections tend to 
zero, in contrast to the ones in Refs. 0, 1^. They also performed a numerical calculation of 
the corresponding Schrodinger equation and compared numerically their formulas with the 
exact values. From the comparison, they have claimed the correctness of their formulas and 
thus the validity of the dilute gas approximation. However, they did not make a comparison 
with our results nor examine the accuracy of them. 

In this letter, considering the present situation described above, we would like to compare 
the formulas by Alhendi and Lashin in Ref. and ours in Ref. 0] from various points of 
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view. In the next section, we first make a numerical comparison of tfie nonperturbative 
corrections. In Section IIIIl we examine the large order behavior of the perturbation series. 
In Section IIVI we take into account the perturbative corrections to check the accuracy of 
the formulas for the ground state. In Section we carry out the ordinary WKB calculation 
for the Schrodinger equation to provide another reference for the comparison. Finally, we 
summarize the results in the last section. 



II. PURELY NONPERTURBATIVE CORRECTIONS 

In Ref. 0, Alhendi and Lashin investigated the following triple- well potential: 

^^^(-;-) = 4^ + Y-V-i)^ (2) 

They calculated the sum of multiple instanton contributions with the dilute gas approxima- 
tion and obtained for the lowest three eigenvalues. 
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On the other hand, the triple-well potential we investigated in Ref. 0] is the following: 

H'^iq; 9^^) = -1^2+ l^'i^ - + ^(1 - SgV). (4) 

Utilizing the valley method, we obtained for all the noth eigenstates localized around the 
central well and the n^-th eigenstates with parity ± localized around the side wells^, 

E^:{9,e) = n. + \ + l + ^.^-''''^ E^}{g^e) + 0(e-/^^), (5a) 

E^Jg, e) = 2n± + 1 - e + 4 ^^'^''^ E^^^g, e) + ©(e^^.^), (5b) 

where the coefficients EnJ and En} are given by 



1 f2\""/ l\~/2-'/4+3./4 . J 3 



^^.S'(9.^) = -rT 3 -4) r(-^ + i-i.), (6a) 



no! Vfi'^y V fl'^y V 2 4 4 

row X ,-l)(i-3^)/2±i /2\'"±+'/'-'^/Vl \"± / 1 3n 



^ Except for the cases e — ±{2J\f + l)/3 {Af = 0, 1, 2, . . .) where a part of the harmonic osciUator spectra of 
the central and side wells degenerates. 
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In order to compare the two results, we must first establish the relation between the Hamil- 
tonians Eqs. © and (^. By applying a scale transformation on the coordinate q in Eq. 
we easily find the following relation: 

H'^^x- u) = ujH''^{uj^/\- 0) . (7) 

Therefore, the quantities which we shall make comparison with Eqs. Q are given by^ 



(8a) 



Ef(uj) = uj^^EZ=^{yj-^'\^) = ^( 1 - ^iJl\-^l''^0[^--) ), (8b) 



Ef{u) = Re i?„^^=o(^"'/', = + ^^/'e'^/^ + 0(e-)^ . (8c) 

Here we note that the nonperturbative corrections in Eqs. (jH)) are in general complex and 
the real parts of them should be taken into account as the spectral shifts; the imaginary 
parts of them are to be canceled with the imaginary parts of the perturbative corrections 
arisingfrom the Borel singularity in the framework of the valley method, see for more details 
Refs. We will later consider the imaginary parts in order to investigate the large order 

behavior of the perturbation series in the next section. 

For the purpose of examining the accuracy of the purely nonperturbative corrections in 
Eqs. dni) and (jHl), it is important to note that in addition to the nonperturbative corrections 
estimated in the formulas there are perturbative corrections to the harmonic oscillator 
spectra^. Hence it does not make sense to compare directly them to the exact eigenvalues 
without taking into account the perturbative contributions. Fortunately, however, the 
perturbative corrections to the harmonic oscillator spectrum for the first and second excited 
states are completely the same (cf. Section HTH) . As a result, the difference of the eigenvalues 
between the first and second excited states IS.E2\ = E2 — Ei only contains the purely 
nonperturbative contributions. Therefore, the comparison of the quantity AE21 enables 
us to study the accuracy of Eqs. Q and (jHI) adequately. In Table HI we show i) the exact 
results Ai?2i obtained by solving numerically the Schrodinger equation for the Hamiltonian 
(121) presented in Ref. ii) ^E^^ = E^^ - Ef"^ obtained from Eq. © and the ratios 
AE^l^/AEIf, and iii) AE^'^ = E^^-Ef^ obtained from Eq. (jH)) and the ratios AE^'^/AEIf. 



From Table I, we see that our results AE2i{uj) are in better agreement with the numerical 
exact values A£'2i(^) than those by Alhendi and Lashin Ai?^'"(co'). We also find that as 
the parameter u becomes larger (the coupling g becomes smaller), the accuracy of our 
results becomes better while that of the results by Alhendi and Lashin becomes worse. This 
suggests that our method provides a reliable semi-classical approximation but the dilute gas 



^ In this case, the spectral sphtting takes place between the first and second excited states and hence the 
parity odd state is lower than the parity even state in the spectrum, which is in contrast to the case of 
symmetric double-well potentials. 

^ In this article, a contribution is called perturbative (nonperturbative) if it can (cannot) be expressed as a 
(formal) power series in g — lu^^^'^, respectively. See also Section Hvl 
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TABLE I: Comparison of the energy difference between the first and second excited states. 
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30 


4.7230 • 10" 


-3 


3.7381 • 


10- 


-3 


0.79147 


6.8061 


10- 


-3 


1.4411 


50 


9.1006 • 10- 


-7 


4.7154 • 


10" 


-7 


0.51814 


1.1081 


10- 


-6 


1.2176 


70 


1.0186 • 10- 


'10 


4.1959 • 


10" 


-11 


0.41193 


1.1667 


10- 


-10 


1.1454 


90 


8.9504 • 10- 


'15 


3.1490 • 


10" 


-15 


0.35183 


9.9282 


10" 


-15 


1.1092 


110 


6.8449 • 10- 


-19 


2.1356 • 


10- 


-19 


0.31201 


7.4439 


10- 


-19 


1.0875 



approximation does not. As is easily seen from the Euclidean action of the present model 

the coupling constant g"^ plays a similar role to the Plank constant h. Therefore, any reliable 
semi-classical approximation should work better as the coupling constant g becomes smaller. 

III. LARGE ORDER BEHAVIOR OF PERTURBATION SERIES 

In spite of the fact that perturbation series are in general divergent and at most asymp- 
totic, they contain much information on the property of the physical quantity under consid- 
eration. An intimate relation between nonperturbative property and large order behavior of 
perturbation series is a typical example In this section, we make an examination from 
this point of view. 

In order to evaluate the perturbative corrections to the harmonic oscillator spectrum of 
the each potential well, it is convenient to begin with the Hamiltonian (@)) and then make 
the transformation indicated by the scale relation ((7j) after the perturbative calculation. 
The coefficients of the perturbation series are systematically calculated with the aid of the 
Bender- Wu method First, the perturbation theory around the harmonic oscillator states 
of the central potential well is set up by decomposing the Hamiltonian into the harmonic 
oscillator part and the remaining part: 

^'^(g; 0) = + - 9'<l' + \9'<f- (10) 

The perturbative corrections to the eigenvalues and eigenfunctions are defined by the fol- 
lowing formal series expansions: 

oo oo oo 

E{g) = J2 ^(^' 9) = e-^^/^ Yl 9'' E (11) 

m=0 k=0 1=0 

where P = (1) for the even (odd) parity states, respectively. For the lowest state, c'^l = 1/2 
and a^2i — alH > 0. Requiring that they satisfy the Schrodinger equation, we obtain a 



X = gq, 



(9) 
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recursion relation for Ci^2ilp and c'^™': 

(4/ + 1 + 2P)aflp - (2/ + 2 + P)(2/ + 1 + P)ag^^,)^p 

k 

_9^[2(fc-l)] ,\2{k~2)\ _ c)Sr- \2m]\2{k-m)] 

^"2(/-2)+P "T "2(«-3)+P ~ ^ 2-^ "2/+P • l-^^J 

Second, the perturbation theory around the harmonic oscillator states of the side potential 
wells is defined by shifting the origin of the coordinate to one of the minimum of the side 
potentials q q ± 1/ g and then decomposing the Hamiltonian: 

H^^iq ± 1/g; 9,0) = -^^^ + ^q' ± 6gq' + ^gV ± + ^gV. (13) 

The perturbative corrections are similarly introduced by^ 

OO OO CO 

m=0 k=0 1=0 

The recursion relation for c'™' and ap' in this case is then given by 
(4/ + 2)afl-(/ + 2)(/ + l)a|5,±12af„3^1 

[k/2] 

+ 13af_-/] ± Qal-'^ + af.^'l = 2 c^'""^ af-'l (15) 

m=0 

For the lowest state, ct"' = 1 and a|°' = for all I > 0. 

On the other hand, as we have mentioned previously, the imaginary parts of the non- 
perturbative contributions are to be canceled with those of the perturbative ones in the 
framework of the valley method. This leads to the following dispersion relation 0, |^ : 

Jo 

This relation enables us to predict the large order behavior of the perturbation series for 
the eigenvalues. For the lowest three states (no = n± = 0) in the present case (e = 0), we 
obtain from Eqs. (0)-® and (fT^ 



Here we note that all the perturbative coefficients of odd powers in g for the spectrum vanish due to the 
parity symmetry of the original Hamiltonian p(J|) . 
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TABLE II: The ratios of the exact values of the perturbative coefficients c^^™] to the predicted 
asymptotic values c^^™!. 



m 


[2m] , [2m] 
^0 / 


[2m] , [2m] 


20 


8946472445 


7797002850 


40 


9493285320 


8904365552 


60 


9665686152 


9268736279 


80 


9750492671 


9451085558 


100 


9800964967 


9560611732 

\j a ly '.y \j vy vy -i a- i ■ y ^ 


120 


0.9834448543 


0.9633690036 


140 


0.9858286759 


0.9685922139 


160 


0.9876123158 


0.9725115700 


180 


0.9889971178 


0.9755611660 


200 


0.9901034160 


0.9780016290 


220 


0.9910075563 


0.9799989041 


240 


0.9917603143 


0.9816636733 


260 


0.9923967744 


0.9830725934 


280 


0.9929419564 


0.9842804381 


300 


0.9934141831 


0.9853273870 



Therefore, we can check the vaUdity of the results (0)^® by comparing the predicted asymp- 
totic forms 

^[2m] 

Eqs. flT7|) with the exact perturbative coefficients cl^™' calculated using 
the recursion relations (fT^ and (fT3j) . 

In Table HTl we show the ratios c'^^^/c'^^^ up to the order m = 300. We easily see that 
the exact values indeed tend to the predicted asymptotic values for both the ground and 
excited states and thus confirm the correctness of our formulas, at least, for their imaginary 
parts. 

IV. INTERPLAY BETWEEN NONPERTURBATIVE AND PERTURBATIVE 
CORRECTIONS 

In Section |nj we have investigated the nonperturbative corrections for the excited states 
and confirmed that our formulas obtained from the valley method calculation are more 
accurate than the ones obtained from the instanton calculation with the dilute gas approxi- 
mation. Although the analysis of the large order behavior in the previous section ensures the 
correctness of the imaginary parts of our formulas for both the ground and excited states, 
it does not necessarily mean that the real part of our formula for the ground state is also 
correct. In order to check the accuracy of the nonperturbative spectral shift for the ground 
state we must resort to other means. 

As we have mentioned previously, there exist perturbative corrections to the harmonic 
oscillator spectrum in addition to nonperturbative ones. Hence, each of the spectrum is 
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formally expressed as 



E{g) = E'~'^ +EM9) + EA9), (18) 

where E^^"* denotes the harmonic oscillator spectrum when g = 0, E^p the purely nonper- 
turbative part which cannot be represented by a formal power series in = uj~^, and Ep 
the remaining perturbative part. The decomposition ()18|) suggests that we can check the 
formulas for the nonperturbative corrections by examining the prediction on the perturbative 
corrections instead. That is, we can regard each of the following quantity as the prediction 
of the each formula on the perturbative contribution to the ground state: 

E^'^iuj) ^ E^,^{uj) - E^'^iu;), E^^{u;) ^ E^,^{uj) - E^^{uj), (19) 

where Eq^{u}) is the exact eigenvalue for the ground state, Eq^(uj) and Eq^{u) are respec- 
tively given by Eqs. (jHaj) and (jHaj) . both of which consist of the harmonic oscillator eigenvalue 
and the predicted nonperturbative contribution. 

TABLE III: Predicted perturbative contributions to the ground state energy. 



UJ 






30 


-0.818251- •• 


-0.8 21854- •• 


50 


-0.7883965537••• 


-0.78839 70101--- 


70 


-0.7763341456 10396 • • • 


-0.7763341456 51123- •• 


90 


-0.77005483676111 02506 • • • 


-0.77005483676111 33127 • • • 


110 


-0.76619756312233715 58989 • • • 


-0.76619756312233715 61068 •• • 



In Table IIIH we show the numerical values of the predicted perturbative contributions 
to the ground state energy for each value of u calculated using Eq. (fTUI) . For the exact 
eigenvalues Eq^{u), we have used again the numerical results shown in Ref. 0]. 

To examine the accuracy of these predictions, we shall evaluate the exact perturbative 
contribution Ep by using the perturbation series. Although the perturbation series is gen- 
erally divergent, as we have already observed in the previous section (cf., Eqs. ()17|) ). the 
asymptotic property of the perturbation series nevertheless ensures that for a sufficiently 
small value of the expansion parameter a partial sum of the first finite terms in the pertur- 
bation series gives an asymptotic value of the perturbative correction^. As a consequence, 
however small the value of the expansion parameter is, there exists a critical order at 
which the absolute value of the perturbative correction |5'^"^c|f'"' | becomes minimum. It is 
apparent that the asymptotic property is lost and replaced by the divergent one when the 
order of the perturbation exceeds the critical order nic. Therefore, the asymptotic values of 



^ We note that the asymptotic property of the perturbation series in general has nothing to do with the 
Borel summabihty. We also note that in the quantum mechanical systems there exist no renormalon 
singularities which originate from the IR and UV divergences in the higher-dimensional quantum field 
theories. 
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the exact perturbative corrections we should read from the perturbation series are given by 
the first finite partial sums up to at most the critical order rric. 

In Table IIV[ we illustrate the numerical values obtained from the first finite partial sums 
of the perturbation series for the ground state. Here we note that from the scaling relation 
(|7j) the perturbative quantity Ep^u) we should take for the Hamiltonian H^^{x;u!) reads 

M M 

EA9) = J29""'^o^^ ^ EpM = J].;^-cri (20) 

m=l m=l 

In Table HVl we show the partial sums up to the order M with rric — 10 < M < rric. The 
critical order rric for all the cases u = 30, 50, 70, 90, and 110 are given by rric = uj/2. 

TABLE IV: The first finite partial sums of the perturbative corrections. 





UJ = 


= 30 




cj = 50 




a; = 70 


M 


2-^171= 


, ,l~m„[2m] 

1 Cq 


M 


S^M i_m [2m] 
Z^m=l ^ "-O 


M 


Z^m=l ^ "-0 


6 


-O.J 


? 21307- • • 


16 


-0.78839 69801- • • 


26 


-0.7763341456 49363- • • 


7 


-O.J 


^ 21522- • • 


17 


-0.78839 69881- • • 


27 


-0.7763341456 49752- • • 


8 


-O.J 


^21641--- 


18 


-0.78839 69939- • • 


28 


-0.7763341456 50061--- 


9 


-O.J 


^21714- • • 


19 


-0.78839 69983- •• 


29 


-0.7763341456 50315- • • 


10 


-O.J 


^21764- • • 


20 


-0.78839 70018- •• 


30 


-0.7763341456 50532- • • 


11 


-O.J 


^21800- • • 


21 


-0.78839 70047 • • • 


31 


-0.7763341456 50723- • • 


12 


-O.J 


^21830- • • 


22 


-0.78839 70072 • • • 


32 


-0.7763341456 50897- • • 


13 


-O.J 


^21856- • • 


23 


-0.78839 70095 • • • 


33 


-0.7763341456 51060- • • 


14 


-O.J 


?21880--- 


24 


-0.78839 70118- •• 


34 


-0.7763341456 51217- • • 


rric 


-O.J 


? 21903- •• 


mc 


-0.78839 70140- •• 


rric 


-0.7763341456 51373- •• 





to = 90 




UJ = 110 


M 


V^Af ,l-m^[2m] 
Z^m=l ^ '-0 


M 


Z^m=l ^ 


36 


-0.77005483676111 32120 • • • 


46 


-0.76619756312233715 61012 • • • 


37 


-0.77005483676111 32318 • • • 


47 


-0.76619756312233715 61022 • • • 


38 


-0.77005483676111 32485 • • • 


48 


-0.76619756312233715 61031 • • • 


39 


-0.77005483676111 32628 • • • 


49 


-0.76619756312233715 61039 • • • 


40 


-0.77005483676111 32755 • • • 


50 


-0.76619756312233715 61046 • • • 


41 


-0.77005483676111 32871 • • • 


51 


-0.76619756312233715 61053 • • • 


42 


-0.77005483676111 32978 • • • 


52 


-0.76619756312233715 61059 • • • 


43 


-0.77005483676111 33080 • • • 


53 


-0.76619756312233715 61065 • • • 


44 


-0.77005483676111 33179 • • • 


54 


-0.76619756312233715 61071 • • • 


rric 


-0.77005483676111 33277 • • • 


rric 


-0.76619756312233715 61077 • • • 



Comparing the results in Table IIVI with the ones in Table IIIIl we easily see that the 
asymptotic values of the perturbative corrections for all the cases are in good agreement 
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with the values Ep^^u) predicted by our formula (|Hajl but apparently deviate from the 
values Ep^{uj) predicted by the Alhendi and Lashin's formula (|Haj). 



We note that the differences between Ep^{uj) and 



M 

m=l 



UJ 



Cq are minimum around 
M ~ rric — 1. It indicates that IieEp{g) in the case has the perturbation series as strong 
asymptotic series [l^, that is, there exist positive real constants C and a so that 



M 



Rei?p(^)- J^Z-c^] 



m=l 



2|M+1 



(21) 



for all M and all g^ > near the origin. Indeed, if it is the case, we have with the aid of 
Eq. (ITTal) and a = 2, 



M 



ReEp{g)-J29 



2mJ_2m] 
-0 



m=l 



< C\M + 7/4)^/^ 



2(M+l)„[2(Af+l)] 



(22) 



where C = C7iT{3/4:)/2^/^. The right hand side is minimum around M ~ rric — 1 by the 
definition of the critical order rric, and thus the above fact would be naturally understood. 

Finally, we would like to mention about the fact that in the parameter region we have 
examined, 30 < a; < 110 or 0.095 ^ g 0.18, the following relation is roughly satisfied: 



mm 



ReENp(^) X 10-\ 



(23) 



Interestingly, we can show that a similar relation is observed generically as far as the sys- 
tem under consideration has a nonperturbative effect. Suppose the following conditions 
are satisfied for smaller values of the coupling constant g"^ involved in the system under 
consideration: 

ImEMg) ~ Cg-^^-'+^-'/'a", ReEMg) = AlmE^g), (24) 

where A, C, and b > are real constants. Then, we can prove the following intriguing 
relation for smaller g^: 



mm 



^2m^[2m] 



26e 
1% 



gReE^p{g) 



For the proof, we first note that the first condition in Eq. implies 

TT 

for larger m. Next, we define a function / by 



^2m^[2m] 



(25) 



(26) 



(27) 



It is evident that for larger integer m the function /(m; g) well approximates the magnitude 
of the mth-order perturbative correction. The derivative with respect to fi reads. 



d_ 



f{fi;g) = f{fi;g)[\n{bg^)+^{f, + iy + l 



(28) 
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where denotes the digamma function. Hence /(/i; g) takes minimum value at fi = fi, n 
satisfying 

\n{bg^)+^{jj + iy + l) =0. (29) 

For smaller value of (7^ <^ 1, we see p, becomes larger. Thus, applying the asymptotic 
expansion of the digamma function : 



we obtain 



V;(z)~lnz-l + 0(^-^), (30) 



bg' - exp ( ^ ^ [1 + Oir')] ) . (31) 

H + u + l v2(/i + z/ + 1) -'/ 

With the aid of the Stirling formula and Eq. ()31|) , we have 

r(/^ + 1/ + 1) ~ ^ei/2(^^2)-(A+.+i/2)g-i/6,2 _^ Oir')] . (32) 

Therefore, the minimum value of /(/i, g), which would provide a good approximation to the 
minimum magnitude of the perturbative correction at the critical order rric ~ fl, is estimated 

as, 

fm-g) = ^-^\c\ (bgyrm + u+i) 

TT 



71 



2he 

A-'gReEMg) [l + Oir')], (33) 



and thus we obtain the relation (j25|) . In our case A = 1 and 6 = 2, and thus the relation 
(I23|) for (7 ~ 0.1 follows. As a consequence, we also find that the next-order nonperturbative 
corrections of order ©((^^Rei^Np) becomes negligible in comparison with the perturbative 
correction for sufficiently small g; from the relation we readily obtain 

O (/Re^Np) ~ O (^gmin j^^^ct^^lj . (34) 

Therefore, the next-order nonperturbative corrections do not affect the analysis for (7 ~ 0.1 
in this section. 



V. WKB CALCULATION 

So far, we have checked the accuracy of the semi-classical calculations of the path-integral 
by comparing with the exact values calculated numerically. In this section, we make a 
comparison in a different way. To this end, we employ another nonperturbative approach 
to derive formulas for the same physical quantities. The method we shall use here is the 
WKB approximation for the Schrodinger equation. In the following, we shall always consider 
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the leading terms of the power expansion in g since we are interested in the quantization 
condition for the nonperturbative contribution. 

Let us consider the more general Hamiltonian Q for all eg^ <^ 1. The system has parity 
symmetry and thus it is sufficient to study the connection condition of the WKB wave 
function only on the half-line q G (0, oo). In the vicinity of the minimum of the central 
potential well, the Schrodinger equation for the Hamiltonian (0]) in the leading order of g is 
given by 

1 d^ 



(35) 



The local solutions possessing a definite parity ± are expressed as 

iP{q) = A± (D,i-V2q) ± D,iV2q) 



(36) 

where A± are constants and is the parabolic cylinder function with u = E — e/2 — 1/2. In 
a similar way, around the minimum of the right side potential well, the Schrodinger equation 
is approximated by 

'2dg2 
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^{q) = {E + e)i,{q). 



(37) 



and the local solution which vanishes at g ^ oo is given by 

,lj{q) = BDx{2{q-l/g)), 
where i? is a constant and A = E /2 + t/2 — 1/2. The solutions ()36|1 and 



(38) 



are to be 

connected with the following WKB solution in the classically forbidden region (gi <^ g -C ^2)- 



il){q) 



A;(g)i/2 



exp 



dx k{x) \ + 



91 



Co 



k{qy/' 



exp 



dx k{x] 



91 



where 



k{x) = v/x2(l - g^x^y + e(l - 3^23,2) _ 2E . 



(39) 



(40) 



The positive classical turning points qi{i = 1,2) with < gi < g2 defined by the solutions 
of V{qi) = E are. 



V2E-e + 0{g'^ 



1 

q2 = - 
9 



E + e 



+ 0{g). 



(41) 



In order to connect the wave functions obtained in the each region, it is important to note 
that the leading term in g of the WKB solution (jH^ll varies according to the position it is 
viewed from. If it is viewed from the point around the central potential well, the integral in 
the exponent in Eq. (jHTIjl is evaluated as 

rq 



dx k{x) 
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\ r du^w\l - w^f + eg\l - 3^2) _ 2E 

1 J C|q^ 

1 {2E - e)g^ ?>eg^ 



9' 



91 



du) 



ggi 



2? 



4^ 



w{l — w ) 

e - 2E 
+ — ^ — 



w 



2 w(l 



2 1 



€ -\- E 

In \iu\ H In |1 — tu^l + ■ 



+ 



1 91 



591 



(42) 
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and k{q) ~ g + ■ ■ ■ . Thus, in the leading order of g we obtain the WKB wave function as 



ip{q) 



,1/2 



{2E-e)/2 ^ 



1/2 



f ^ 



— e 



-(2£;-e)/2 



(43) 



Comparing this with the following asymptotic form for g 3> 1 of the wave function ()36|) (cf. 
Ref. ^Ij) determined in the region of the central potential well: 



Hi) ~ A 



± 



2n 



{i-ir ± l)e-^/^(v^g)^ + -^e'V^(v^g)— ^ 



we have the following connection condition: 

\/2^ _C2 ( 2y/E-€/2' 



(44) 



2E- 



{{-ir±i)T{-u) c. 



(45) 



On the other hand, the integral in the exponent in Eq. ()39|1 viewed from the point around 
the right side potential well is evaluated as 



dx k{x) 



~2 



91 











1 








~7 



9 Jgqi-l 

duj 



, eg\2 + 6uj + 3u;^) + 2Eg^ 



2uj{l + uj){2 + uj) 



UJ 

T 



+ 



e + E 



ln|a;(2 + u;)| 



e-2E 



lull + Lc;| + ■ ■ 



(46) 



and k{q) ~ 1/ g — q + ■ ■ ■ . Thus, in the leading order of g we obtain the WKB wave function 
as 



+ 



{1/g-qy/' 

C^Qy^9'g3E/2 



{2E-t)/2 



<^-'l^)\2{l/g-q)) 



\V2E^e, 

/ N -(2S-e)/2 

\V2E^e. 



(47) 



Matching this with the following asymptotic form for 1/g — q ^ 1 of the wave function 
determined in the region of the right side potential well: 



^(g) ~ B 



'(-V^)^(2(g - 1/,))^ - ^^^^^^e(-V.)^(2(g - 1/,)) 



A-l 



r(-A) 



(48) 



we get another connection condition as follows 

V2^ 



-i)^+ir(-A) C2 



V2E 



2E-e 



(49) 
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Therefore, eliminating the coefficient C1/C2 in Eqs. fl45p and ()49|) . we finally obtain the 
following quantization condition: 




Strikingly, this condition is in complete agreement with the one obtained previously by us 
with the valley method, Eq. (5.19) in Ref. from which our formulas of the nonperturbative 
effect (0)-® have been derived. We note that the coincidence is achieved not only for the 
case e = we have examined in the preceding sections but also for all e 7^ where we cannot 
apply the usual instanton technique since the classical configuration now becomes a bounce 
solution. 



VI. SUMMARY 

In this letter, we have extensively examined lowest three energies of the symmetric triple- 
well potential with non-equivalent vacua by comparing the two different sets of the formulas, 
the one is calculated by means of an instanton technique with the dilute gas approximation 
in Ref. 0] and the other is by the valley method beyond this approximation scheme in 
Ref. 0. 

First, we compared numerically both the formulas with the exact values for the spectral 
splitting between the first and second excited states due to the purely nonperturbative effect. 
We have found that in contrast to the latter formulas the prediction of the former formulas 
get worse as the value of the parameter tends to the region where the semi-classical methods 
should work better. Thus, contrary to folklore that this kind of problem can be handled 
by the use of the dilute gas approximation, our careful comparison with the exact results 
indicates that the dilute gas approximation is insufficient to produce the correct asymptotic 
formulas even if we restrict ourselves to examining the low lying eigenvalues. 

Furthermore, we examined for the ground state both the perturbative and nonperturba- 
tive contributions. We have found that the asymptotic values of the perturbative corrections 
calculated from the perturbation series deviate from the predicted values obtained from the 
instanton calculation with the dilute gas approximation but are in good agreement with the 
ones obtained from the valley method calculation. 

We also checked the accuracy of the imaginary parts of the latter formulas by testing the 
prediction on the large order behavior of the perturbation series. We evaluated the pertur- 
bative coefficients up to the order 300 and compared them with the predicted asymptotic 
behaviors. We have confirmed the correctness of the prediction for both the ground and 
excited states. 

Finally, we carried out the WKB calculation in the leading order of the coupling constant. 
The resulting WKB quantization condition is in complete agreement with the one obtained 
by means of the valley method in Ref. from which our formulas have been derived. In 
other words, the dilute gas approximation in the path integral formalism does not correspond 
to a proper semi-classical approximation. In this respect, we would like to recall the fact 
that discrepancy between the dilute gas approximation and the WKB approximation has 
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been already reported in Ref. 3| for the wave functions even in the case of the symmetric 
double-well potential. 

To conclude, all the present analyses entirely support the valley method calculation and 
indicate the limitation of the dilute gas approximation in the present triple-well potential 
problem. Therefore, we would like to repeat the assertion in Ref. ^ that the applicability 
of the dilute gas approximation would be quite hmited. 
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